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EXACT  SOLUTIONS  OF  MODERATELY  THICK  LAMINATED  SHELLS 
By  J.  N.  Reddy,1  M.  ASCE 


ABSTRACT: 

The  paper  contains  an  extension  of  Sanders  shell  theory  for  doubly-curved 
shells  to  a  shear  deformation  theory  of  laminated  shells.  The  theory  accounts 
for  transverse  shear  strains  and  rotation  about  the  normal  to  the  shell  mid¬ 
surface.  Exact  solutions  of  the  equations  are  presented  for  simply  supported, 
doubly-curved ,  cross-ply  laminated  shells  under  sinusoidal,  uniformly  distrib¬ 
uted,  and  concentrated  point  load  at  the  center.  Fundamental  frequencies  of 
cross-ply  laminated  shells  are  also  presented.  The  exact  solutions  presented 
herein  for  laminated  composite  shells  should  serve  as  bench  mark  solutions  for 
future  comparisons. 

INTRODUCTION 

There  exist  a  number  of  theories  for  layered  anisotropic  shells.  Many  of 
these  theories  were  developed  originally  for  thin  shells,  and  are  based  on  the 
Kirchhoff-Love  kinematic  hypothesis  that  plane  sections  normal  to  the  unde¬ 
formed  midsurface  remain  plane  and  normal  to  the  middle  surface  after  deforma¬ 
tion  and  undergo  no  thickness  stretching.  Surveys  of  various  shell  theories 
can  be  found  in  the  works  of  Naghdi  [161  and  Bert  [3,41,  and  a  detailed  study 
of  thin  ordinary  (i.e.,  not  laminated)  shells  can  be  found  in  the  monographs 
by  Kraus  [141,  Ambartsumyan  Fll,  and  Vlasov  [201. 
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The  first  analysis  that  incorporated  the  bending-stretching  coupling  (due 
to  unsymmetric  lamination  in  composites)  is  due  to  Ambartsumyan  [2].  In  his 
analyses  Ambartsumyan  assumed  that  the  individual  orthotropic  layers  were  ori¬ 
ented  such  that  the  principal  axes  of  material  symmetry  coincided  with  the 
principal  coordinates  of  the  shell  reference  surface.  Thus,  Ambartsumyan's 
work  dealt  with  what  is  now  known  as  laminated  orthotropic  shells  rather  than 
laminated  anisotropic  shells;  in  laminated  anisotropic  shells  the  individual 
layers  are,  in  general ,  anisotropic  and  the  principal  axes  of  material  symme¬ 
try  of  the  individual  layers  coincide  with  only  one  of  the  principal  coordi¬ 
nates  of  the  shell  (the  thickness  normal  coordinate). 

Dong,  Pister,  and  Taylor  [9l  formulated  a  theory  of  thin  shells  laminated 
of  anisotropic  material  that  is  an  extension  of  the  theory  developed  by 
Stavsky  [91  for  laminated  anisotropic  plates  to  Donnell's  shallow  shell 
theory  (see  Donnell  [10]).  Cheng  and  Ho  [8]  presented  an  analysis  of  lamina¬ 
ted  anisotropic  cylindrical  shells  using  Flugge's  shell  theory  (see  Flu'gge 
[121).  A  first  approximation  theory  for  the  unsymmetric  deformation  of  non- 
homogeneous,  anisotropic,  elastic  cylindrical  shells  was  derived  by  Widera  and 
his  colleagues  [25,26]  by  means  of  the  asymptotic  integration  of  the  elastici¬ 
ty  equations.  For  a  homogeneous,  isotropic  material,  the  theory  reduces  to 
Donnell 's  equations. 

All  of  the  works  reviewed  above  are  based  on  Kirchhoff-Love's  hypotheses 
in  which  the  transverse  shear  deformation  is  neglected.  These  theories,  known 
as  the  Love's  first-approximation  theories  (see  Love  [151)  are  expected  to 
yield  sufficiently  accurate  results  when  (i)  the  lateral  dimension-to-thick- 
ness  ratio  is  large;  (ii)  the  dynamic  excitations  are  within  the  low-frequency 
range;  and  (iii)  the  material  anisotropy  is  not  severe.  However,  application 
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of  such  theories  to  layered  anisotropic  composite  shells  could  lead  to  as  much 
as  30*  or  more  errors  in  deflections,  stresses,  and  frequencies. 

The  effects  of  transverse  shear  deformation  and  transverse  isotropy,  as 
well  as  thermal  expansion  through  the  thickness  of  cylindrical  shells  were 
considered  by  fiulati  and  Essenberg  f 121  and  Zukas  and  Vinson  [291.  Whitney 
and  Sun  [23,241  developed  a  shear  deformation  theory  for  laminated  cylindrical 
shells  that  includes  both  transverse  shear  deformation  and  transverse  normal 
strain  as  well  as  expansional  strains. 

Recently,  Bert  and  his  colleagues  [6,7,13]  presented  exact  solutions  for 
bending  and  vibration  of  two-layer,  cross-ply,  thin  cylindrical  shells.  These 
solutions  are  limited  to  cylindrical  shells  and  sinusoidal  distribution  of  the 
transverse  load,  and  the  procedure  used  is  one  similar  to  that  used  by  Whitney 
and  Leissa  [21],  Whitney  and  Pagano  [22],  Bert  and  Chen  [5],  and  Reddy  and 
Chao  [17]  for  laminated  composite  plates.  The  present  study  is  concerned  with 
the  development  of  exact  solutions  for  simply  supported,  doubly-curved,  cross- 
ply  laminated  shells.  It  is  shown  that,  unlike  plates,  antisymmetric  angle- 
ply  laminated  shells  with  simply  supported  boundary  conditions  do  not  admit 
exact  solutions.  Numerical  results  for  transverse  deflection  and  fundamental 
frequencies  are  presented  for  various  cross-ply  laminated  shells. 


DERIVATION  OF  EQUATIONS 

Figure  la  contains  a  differential  element  of  a  doubly-curved  shell. 

Here  U|,52,C)  denote  the  orthogonal  curvilinear  coordinates  (shell  coordi¬ 
nates)  such  that  the  £j-  and  fy-  curves  are  lines  of  curvature  on  the  midsur¬ 
face  C*0,  and  C-  curves  are  straight  lines  perpendicular  to  the  surface  c=0. 
For  the  doubly-curved  shells  discussed  here,  the  lines  of  principal  curvature 
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coincide  with  the  coordinate  lines.  The  values  of  the  principal  radii  of 
curvature  of  the  middle  surface  are  denoted  by  Rj  and  R£. 

The  position  vector  of  a  point  on  the  middle  surface  is  denoted  by  r,  and 
the  position  of  a  point  at  distance  z  from  the  middle  surface  is  denoted  by  R 
(see  Fig.  lb).  The  distance  ds  between  points  and  (^+d£ , 

is  determined  by 

(ds)2  =  dr  •  dr 


=  a|(d^)  2+  a^(dZp)2  (1) 

Sr 

where  dr  =  r^d^  +  r^d^ »  the  vectors  and  r\>  (r^  =  are  tangent  to 
the  Si  and  zp  coordinate  lines,  and  and  a-,  are  the  surface  metrics 


<*?  =  JCi  *  Ci  .  <*$  s  Up  *  C2* 

The  distance  dS  between  points  and  ( Si  +d  ,  Ip >  C+d  C )  is  given  by 

(dS)2  =  dR  •  dR 


=  L^dSi)2  +  L  ^(dS?)2  +  L|(dc)  2 


(3) 


dR 


dR 


dR 


where  dR  =  ds^  +  d^  +  dc,  and  Lj,  L£,  and  L3  are  the  Lame1 
coefficients 

4  =  “l  ^  +  R^  ’  4  =  ^  +  *  L3  =  1 


(4) 


dR  dR 

It  should  be  noted  that  the  vectors  and  are  parallel  to  the  vectors 
dr  dr  ar’2 


— —  =  r  and  — -  =  r0. 
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From  Fig.  la,  the  elements  of  area  of  the  cross  sections  are 


dA1  =  Ljdqdc  =  otjfl  + 


dAg  =  L  2^  ^2^ ^  =  ^2^  +  R 


(5) 


Let  be  the  tensile  force,  measured  per  unit  length  along  a  £<>-  coor¬ 
dinate  line,  on  a  cross  section  perpendicul ar  to  a  coordinate  line.  Then 
the  total  tensile  force  on  the  differential  element  in  the  ^-direction  is 
N^apd^.  This  force  is  equal  to  the  integral  of  o^dA^  over  the  thickness, 

h/2 

Ni  a^d^p  =  /  <j|dAj>d^  (6) 


where  h  is  the  thickness  of  the  shell  (C  *  -h/2  and  C  =  h/2  denote  the  bottom 
and  top  surfaces  of  the  shell).  Using  Eq.  (S),  we  can  write 


h/2  r 

N1  =  /  »j(l  +  R  )dC 
-h/2  2 


(7) 


Similarly,  the  remaining  stress  resultants  per  unit  length  can  be  derived. 
The  complete  set  is  given  by  (see  Fig.  1c). 


Ni  ' 

f  ojd  +  c/«2)  ' 

n2 

0p(l  +  C/R^ 

*12 

ofi(l  +  C/R?) 

"21 

h/2 

o6(l  +  C/Rj) 

N, 

■  Oi 

r  f 

-h/2 

o5d  +  c/r2) 

02 

a4(l  +  C/Rj) 

■******~-^^^^-- - - r 
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Mi 

Co1(l  +  C/R?) 

m2 

Co2(l  +  C/R ! ) 

m12 

Co6(l  +  C/R2) 

m21 

N 

Note  that,  in  contrast  to  the  plate  theory  (which  is  obtained  by  setting 
1/Rj  =  l/R^  =  0),  the  shear  stress  resultants  N^,  and  M21,  and  twisting  mo¬ 
ments  M12  and  M21  are,  in  general,  not  equal.  For  shal 1 ow  she! 1 s ,  however, 
one  can  neglect  c/R^  and  C/R2  comparison  with  unity.  Under  this  assump¬ 
tion,  one  has  and  N?1  =  ar)(i  ^12  anci  M21  s  course»  for  spherical 

shells  we  always  have  Nj?  =  N2^  and  M^2  =  M2^). 

A  set  of  simplifying  assumptions  that  provides  a  reasonable  description 
of  the  behavior  of  thin  elastic  shells  is  used  to  derive  the  equilibrium  equa¬ 
tions  that  are  consistent,  via  virtual  work,  with  the  assumed  displacement 
field: 

1.  the  thickness  of  the  shell  is  small  compared  to  the  principal  radii 
of  curvature  (h/R^  h/R2  «  i); 

2.  the  transverse  normal  stress  is  negligible; 

3.  normals  to  the  reference  surface  of  the  shell  before  deformation 
remain  straight  but  not  necessarily  normal  after  deformation  (a 
relaxed  Ki rchhoff-Love ' s  hypothesis); 

4.  The  deflections  are  small  compared  with  the  shell  thickness. 

The  shell  under  consideration  is  composed  of  a  finite  number  of  ortho¬ 
tropic  layers  of  uniform  thickness,  as  shown  in  Fig.  2.  In  view  of  assumption 
l,  the  stress  resultants  in  Fq.  (8)  can  he  expressed  as 


% 
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N  ck 

(N-i  )  =  I  /  a.(l,C)dC  ,  1  •  1,2,6 

k=1  <*-1 

N  2  Ck 

Oj  s  E  Kf  I  a;dC  ,  i  =  4,5  (9) 

k  =  l  1  r 

1  ck-l 

where  N  is  the  number  of  layers  in  the  shell,  and  are  the  top  and 

bottom  c-coordinates  of  the  k-th  lamina,  and  are  the  shear  correction  factors 
The  strain-displacement  equations  of  a  shell  are  an  approximation,  within 
the  assumptions  made  above,  of  the  strain-displacement  relations  referred  to 
orthogonal  curvilinear  coordinates.  In  addition,  we  assume  that  the  trans¬ 
verse  displacement  u3  does  not  vary  with  C.  r  in  the  shear  deformable  theory 
of  flat  plates,  we  begin  with  the  displacement  <eld 

“1  =  (Llul^  +  C*1  ’  "2  =  ^  (4u2^  +  ^  ’  “3  =  u3  (ir» 

where  (u^.u^.u^)  are  the  displacements  of  a  point  (^,S2,C)  along  the 
Uj ,S?,C)coordinates  and  (ui,u2,u3)  are  the  displacements  of  a  point 
(£j,);p,0).  A  higher-order  (in  C)  displacement  field  can  he  assumed  in  place 
of  Eq.  (10).  Such  theories  have  been  considered  by  Uhitney  and  Sun  [231.  In 
the  interest  of  brevity,  although  the  procedures  to  he  described  can  he  ap¬ 
plied  to  any  higher-order  theory,  we  consider  only  this  first-order  shear  def¬ 
ormation  theory.  Substituting  Eq.  (10)  into  the  strain-displacement  relations 
of  a  orthogonal  curvilinear  coordinate  system,  and  assuming  that  ai  are  con¬ 
stant  (to  specialize  the  results  to  doubly-curve  shells)  one  obtains 

el  =  £?  +  <>1 

*2  s  4  +  c<? 
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where 


e4  “  e4 
£5  '  £S 
£6  s  ££  +  c<6 


1  3ui 

,  o  =  i - 1  +  _i 

1  a1 


i  3u  ■j  u  -j 

-0  =  i _ L  +  _i 

2  cc,  ^ 


n  1  dU3  _  ^2 

^  <*2  ^  R? 


i  ou  o  u  i 

e5  =  ^  +  *1  '  ^ 


n  _  1  5u?  (  1  5ul 

'6  acj  d^2  a?  d^p 


1  3*1 


<1  al  5"l 


K2  d?  3C; 


al  ^^1  '*7.  '^‘’2 


1  9<tl  ,1/1  1  wl  &U? 

Oj  5^  9  ^2  a2 


1_  ?U 1 ^ 

flt^3  5 


where  and  t>9  are  the  rotations  of  the  reference  surface  ^=0  about  the  T0  - 

and  ^  -  coordinate  axes,  respectively.  It  should  be  noted  that  the  displace¬ 
ment  field  in  Eq.  (10)  can  be  used  to  derive  the  general  shear  deformation 
theory  of  laminated  shells  (i.e.,  are  not  necessarily  constant). 

The  stress-strain  relations  for  the  k-th  ortnotropic  lamina  i"  the  mate¬ 
rial  coordinate  axes  are  given  by 


O'-  >  w'  -l/'  * 
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f  -  1 

— 

— 

al 

C11 

c12 

0 

0 

0 

c12 

£22 

0 

0 

0 

> 

0 

0 

c44 

0 

0 

0 

0 

0 

c55 

0 

Or 

6 

(k) 

0 

0 

0 

0 

c66 

k  J 

(k) 


£1 
£2 
1  e4 
e5 
£6 


(13) 


(k) 


where  c^j  are  the  (plane  stress-reduced)  material  stiffnesses  of  the  lamina: 


.  c, 


E2  v12 


11  v12  V21  12  ^~v12^21 


,  c 


22  1  *  ^12 ^2 1 


c44  =  G23  »  c55  =  G13  ♦  c66  =  G12 


(14) 


and 


Ej.E;?  =  Young's  moduli  in  1  and  ?  materi  al -principal  directions, 

respectively 

v^j  =  Poisson's  ratio  for  transverse  strain  in  the  j-th  direction 

when  stressed  in  the  i-th  direction 
G?3,G13,G12  =  s^ear  moduli  in  the  2-3,  1-3,  and  1-2  surfaces, 
respectively. 

Poisson's  ratios  and  Young's  moduli  are  related  hy  the  reciprocal  relations 

v-j  jEj  =  vj  ^  E  ^  (i  =  1,2).  (15) 

The  stress-strain  relations  (13)  transformed  to  the  shell  coordinates  become 


M  -  [o(k)]M 


(16) 


,(k) 


where  are  the  material  properties  of  k-th  layer. 
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The  principle  of  virtual  work  in  the  present  case  yields 


h/2 


0  =  /  If  [a^<5e1+4k^5e?+4k'6E6+o^5Eil+05^6E5]a1a;,dr3ld^}dC 

-h/2  Q 


-  /  [N^Se^  +  +  Ng +  M^6ic^  +  ^2^K2  +  ^g<!>Kg  + 


+  C^SeiJ  -  qSu^  jaj^d^d^ 


(17) 


where  q  is  the  distributed  transverse  load. 

The  governing  equations  of  equilibrium  can  be  derived  from  Eq.  (17)  by 

o 

integrating  the  displacement  gradients  in  e^  by  parts  and  setting  the  coeffi¬ 
cients  of  5u .j  (i  =  1,2,3)  and  6 (i  =  1,2)  to  zero  separately.  Thus  one 
obtains 

3N,  f), 

+  _2_  fN  +  C  NL  )  +  -i  =  0 
ax1  dx^  1  6  0  6'  Rj 

^  dNr,  Q9 

—  (N,  -  c  «,)  +  — =  0 
dXj'6  o  51  3x ,,  R ,, 

+  ^2  r!li  + !!?  ,  =  n 

3x^  3X2  -Rj  Ro  d 

3M,  QMr 

— -  +  — ; -  -  0.  =0 
3x1  3x?  1 


3Mft  3M„ 

1*7 +  w2  -  °2 =  0 


(18) 


where  cQ  denotes  the  constant 


'0  2 


I  fi 

o  tR 


1 


R?^  ’  d*i 


“id^i 


(19) 
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This  term  is  introduced  by  Sanders  (see  [18")),  and  distinguishes  the  Sanders 
theory  from  others.  Using  Eq.  (13)  in  Eq.  (8),  the  stress  resultants 
(Ni»Mi,01)  can  be  related  to  ( e°,<i )  by 

Ni  =  AijeJ  +  ®i  jKj 

(i,j  =  1,2,6) 

Mi  =  8ijej  +  °ij,fj 


Q2  =  A44e4  +  A45  e5 
Q1  =  A45e$  +  A55£5> 


(20) 


Here  A^,  and  ( i , j  =  1,2,6)  denote  the  extensional,  flexural- 

extensional  coupling,  and  flexural  stiffnesses: 

N  ^k 

(Aij*Bij‘°ij)  =  l  /  Qij}  (1 *C.C2)dc-  (21) 

k=1  ck-l 

The  equations  of  Love's  first-approximation  shell  theory  can  be  obtained 
by  setting  c0  =  0  in  Eq.  (18).  Equations  (18)  can  be  specialized  to  flat 
plates,  cylindrical  shells,  and  spherical  shells,  respectively,  by  setting 

=o  and  R?  -  R  (the  x,-axis  is  taken  along  a  generator  of  the 

.r<i  Hp  * 

cylinder),  and  Rq  =  R2  =  R*  The  classical  thin  shell  theory  can  be  obtained 
du3  Uj  du3 

by  setting  *j-  -  ,  and 

CLOSED-FORM  SOLUTIONS 

The  exact  form  of  the  spatial  variation  of  the  solution  of  Eqs.  (18)  can 
be  obtained  under  the  following  conditions: 
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(i)  Symmetric  or  antisymmetric  cross-ply  laminates:  i.e.,  laminates 
with 

A16  =  A26  =  B16  =  B26  =  °16  =  °26  =  A45  =  °* 

(ii)  Freely  supported  boundary  conditions: 

Nj(0,X2)  =  (a,X2 )  =  M^(0,X2)  -  M^(a,X2)  =  0 

03(0, X2)  =  U3(a,X2)  =  U2(0,X2)  =  U2(a,X2)  =  0 

N?(xx  ,0)  =  N2<xlfb)  =  MjjtXj.O)  =  ftjfxj.b)  =  0 

u3(xi »0)  =  u3(x1,b)  =  u1(x1,0)  =  Ujfxj.b)  =  0 

=  <f2(a,x  2)  =  4>^(*j»0)  =  ^^(xj.b)  =  0  (23) 

(iii)  Sinusoidal  (spatial)  distribution  of  the  transverse  load: 

CD 

q  •  l  q  sinax.sinpXg  ,  a  =  J2  .  P  =  jp  (24) 

m,n 

where  a  and  b  are  the  dimensions  of  the  shell  middle  surface  along 
the  x^  and  X2  axes,  respectively. 

The  exact  form  of  the  spatial  variation  of  u-j  (i  =  1,2,3)  and  (i  =  1,2)  is 
given  by 

u1(x1,x2)  =  l  Umncosax1sinpx2 
'  m,n 

u2(xi,x2)  =  l  Vninsinax1cospx2 
m,n 

u3 (X1  ,x2  ^  =  £  wmnsinax1sinpX2 

m,n 

^(xj.Xg)  =  T  Xmncosax1sinpX2 
m,n 

( x  1 » x  2^  =  £  Ymnsinax1cospx2.  (25) 

m,n 

Clearly,  the  solution  satisfies  the  boundary  conditions  in  Fq.  (23). 
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Substitution  of  Eqs.  (24)  and  (25)  into  Eq.  (18)  yields  a  set  of  five 
linear  algebraic  equations  in  terms  of  the  unknown  amplitudes  U,  V,  W,  X,  and 
Y.  These  equations  can  be  expressed  in  matrix  form  as 


for  bending,  and 


-  [Cl {A}  =  {F} 


-  CC1  {A}  •  cj2[S  1  {a} 


for  natural  vibration.  Here  to  denotes  the  frequency  of  natural  vibra¬ 
tion,  {a}  and  {F}  denote  the  columns 


{a}T  =  {U  ,V  ,W,X  ,Y }  ,  (F}T  =  {0,0,qmn,0,0}, 


and  c-j  =  Cji  and  =  S j ^  are  given  by 


C11  =  "  “2A11  “  p2A66  '  C02p2°66  '  TT  ’  2cop2R66 

K1 


c12  =  -  <x0(A12  +  A66  -  c  2d66  ) 


.  .  fAll  ,  A12  A55 i 

c13  -  a  R2  \ 


=  -  a2Rn  -  32B6fi  -  CQf?2D^fi  +  R 


c15  *  -  ae(Bl2  +  366  +  con6fi) 


c22  =  'a\6  ‘  a2co°fi6  "  p2A22  _  7 T  +  2coa2B66 

k2 
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where 


being  the  density  of  the  k-th  layer, 


QrA12  .  A22  .  ^44  \ 
"  ^  R1  +  R2  +  R2 


•23 


c24  =  -  ap(R66  -  con66  +  812) 


44 


c25  =  “  °^B66  *  co°66)  "  ^B2?  + 


'33 


A55*2 


A44^2 

1 

'  R1 

Aii 

A12  \  1 
V] '  ‘ 

c34  = 

rBn 

r12 

R2 

"  A55  ^ 

C35  = 

2 

B22 

+  R 

R2 

-  A44  ^ 

44  =  " 

a2°ll 

-  P2», 

66  '  A55 

■c45 

=  -  ap(n12  + 

°6fi) 

,A12  V? 

IVD  0 

'2  K1  2 


c55  =  '  a2°66  '  ®2°22  '  A44 

Su  =  Pj  +  2P2/R1,  S14  =  P2,  S22  =  Px  *  2P2/R2,  S25  =  P? 


S23  =  Pi»  S44  =  S55  =  Pv  all  other  Sij  =  0. 
N  C, 


(P 


1,R?,P3^  *  l  P(k)U»C,C2)dC 
1  k-l  Ck_! 


(29) 


(30) 


(31) 


Similar  calculation  for  antisymmetric  angle-ply  shells  (A16  =  A26  =  Bu 
B12  “  B22  =  8 65  =  °16  =  °26  =  A45  s  °)  with  the  so1ution  of  the  form 


1 

A 
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ul(xl’x2>  =  I  UmnsinaXlcosPx2 
m,n 

u2(xi  ,x2)  =  l  Vmncosax1sinpx? 
m ,  n 


(32) 


and  U3,  4>1 ,  and  ^  as  given  in  Eq.  (25),  shows  that  the  exact  solution 
is  possible  if  and  only  if  1/Ri  «  1/R2  =  0-  I"  other  words,  antisymmetric 
angle-ply  plates  admit  exact  solutions  but  shells  do  not.  To  see  this,  sub 
stitute  Eq.  (24),  last  three  of  Eq.  (25),  and  Eq.  (32)  into  Eq.  (18).  Note 
that  the  assumed  solution  satisfies  the  following  type  of  simply  supported 
boundary  conditions. 

Ns(0»x?)  =  f%(a,X2)  =  Ng(xpO)  =  Ng(x^,b)  =  0 

u^O.Xg)  =  u1(a,x2)  =  u2(x1,0)  =  u2(x1,b)  =  0 

u3(°,xi)  =  u3(a,x2)  =  u3(x1,(l)  =  u3(x1,b)  =  0  (33) 

$l(xl*0)  =  Cx  1  *  ^  )  =  =  <t>2  (a » xj> )  =  3 

M2(x1§0)  *  M2(xj,b)  =  MjfO.Xp)  =  Mj(a,x2)  =  0 


From  the  first  equation  in  Eq.  (18),  we  obtain 


A4, 

{[-2<xP(Alfi  +  C0B16)]U(nn  +  [-a2(A16  -  coBlfi)  -  P2fA?g  +  co8?f0  -  R^-]vmn 

f  K„  *  *  cno66,  *  ^!*„n 

*  [-«e(>l2  *  »66  +  cor,66')<mnlc0S«XiSlnflx2 

*  (t-»2»n-  P2(»66t  Jcoa68*  co°66>  -  *  [-e(A12*  2coS66*  Co°66' lVmn 

A.  f-  +  c  B1(  o_Boc  A-r 

+  [  16  0  —  +  — —  -  2-6-  +  -4i]gW  +  [-ap(2B.fi  +  c  0  )]X 

L  R^  R2  Rj  ^  mn  1  16  0  16  ' J 

A45 

+  [  ~ a2B  ig“  P2  (B^6  +  co°26  ^  +  — ]Ymn^si  naXlC0Spx2  =  0 


nn 


(34) 


IB 


For  antisymmetric  angle-ply  laminates,  the  following  stiffnesses  are 
identically  zero: 

A16  =  A26  3  811  =  S12  =  B22  3  B66  =  D16  3  °26  3  A45  =  0  (35) 


As  a  result  Eq.  (34)  reduces  to 


A11  .  A12  .  A55 ■ 


[-  2a0coB16Umn  +  co(a2B16  -  P2B2fi)Vmn  +  at(q  +  »  +  r~ )w. 


1  2  '1 


+  82co°6e  R^Xmn  "  a^VnJCOSaXlsin8x2 


+  {[-  «2An  -  P2(A66  +  co°66^  ■  r2  JUmn  +  C”  a0(Al?  +  A66  +  coD%)]vn 

B  8 

+  •  2a**16X«n  -  (aJR16  *  ^26  >Vm,,  >  ■ 


Thus,  we  get  two  equations,  by  setting  the  coefficients  of  cosaXjSinfix-,  and 
sinax^cospx^  to  zero,  from  the  first  equilibrium  equation.  Similarly,  we  ob¬ 
tain  eight  more  equations  from  the  remaining  four  equilibrium  equations.  In 

W„ ... 

other  words,  we  have  ten  equations  in  five  unknowns  (|Jmn,  Vmn,  Xmn,  Ymn), 

r- 

which  have  no  unique  solution.  These  ten  equations  reduce  to  five  when  1/Rj 
and  1/R2  are  set  to  zero,  as  can  be  seen  from  Eq.  (36).  Thus  antisymmetric 
angle-ply  shells  do  not  admit  exact  solution  while  antisymmetric  angle-ply 
plates  do. 
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NUMERICAL  RESULTS 

As  a  first  example  of  bending,  a  doubly  curved  shell  under  transverse 

and 

concentrated  load  at  the  center  is  analyzed.  The  following  geometriC|material 
parameters  are  used: 


Rl  =  R2  =  96.0  in.,  a  =  b  =  32.0  in.,  h  =  0.1  in., 

E^  =  Eg  -  107  psi,  v  =  0.3,  intensity  of  load  =  100  lb.  (37) 

This  problem  was  also  solved  using  the  finite  element  method  by  Yang  [271. 
Table  1  contains  the  center  transverse  deflection  obtained  using  the  shear  de- 
formationtheory  (SOT)  and  classical  shell  theory  (CST)  for  various  terms  in 
the  series.  The  numerical  solution  of  Vlasov  [20]  is  taken  from  Yang's  paper. 
It  should  be  pointed  out  that  both  Vlasov  and  Yang  did  not  consider  transverse 
shear  strains.  It  is  clear  from  the  results  that  the  series  solution  converg¬ 
es  very  slowly.  The  difference  between  the  values  predicted  by  SOT  and  CST 

theories  is  not  significant  because  the  ratios  a/h  =  320  and  R/h  =  960  are 

very  large  (hence,  the  shell  is  essentially  very  thin  and  shallow). 

To  investigate  the  effect  of  transverse  shear  strains  on  the  center  de¬ 
flection,  the  same  spherical  shell  problem  with  point  load  at  the  center  or 
uniformly  distributed  loading  are  analyzed,  and  the  results  are  presented  in 
Table  2.  Note  that  for  uniformly  distributed  load,  the  solution  given  by  the 
50-term  series  is  the  same  as  that  given  by  the  250-term  series,  indicating 
that  the  convergence  is  achieved  with  50  (or  less)  terms  in  the  series.  Com¬ 
parison  of  the  SOT  results  with  the  CST  results  show  that  the  shear  deforma¬ 
tion  is  significant  for  side  to  thickness  ratios  smaller  than  10.  For  exam¬ 
ple,  the  classical  theory  solution  differs  from  the  shear  deformation  theory 
solution  by  0.357  for  a/h  =  20,  3.67.  for  a/h  =  19,  and  1A.67  for  a/h  =  5. 


TABLE  1.  Center  transverse  deflection  of  a  spherical  shell  under  point  load 
at  the  center  (R^  =  R2  =  96  in:  a  =  b  =  32  in.,  h  =  0.1, 

E  =  107  psi,  v  =  0.3  ,  load,  P  =  100  lbs) 


Theory 

n*=9 

n=49 

n=99 

n=149  n=199  n=249 

SOT 

CST 

VI asov 
[20] 

0.032594 

0.039469 

0.03972 

0.03959 

0.03956 

0.039786  0.039814  0.039832 

0.039647  0.039653 

★ 

nxn-term  series;  finite-element  solution  of  Yang  [271:  0.03867 


TABLE  2.  Center  transverse  deflection  (w  x  103)  of  spherical  shells  under 

point  load  at  the  center  and  uniformly  distributed  load  (see  Table  1 
for  the  problem  data). 


h 

n=49 

n=99 

n=149 

n=199 

n  =  249 

Point  load  at 

the  center 

0.32 

1.60 

3.20 

6.40 

3.664 

0.1646 

0.0349 

0.0067 

3.9019 

0.1713 

0.0376 

0.0080 

3.9194 

0.1735 

0.0386 

0.0085 

3.9270 

0.1748 

0.0393 

0.0088 

3.9319 

0.1757 

0.0397 

0.0091 

3.9356 

0.1764 

0.0400 

0.0092 

Uniformly  distributed  load 

0.32 

314.28  # 
(314.33) 

313.87 

(313.93) 

313.86 

(313.93) 

313.86 

(313.93) 

313.86 

313.86 

1.60 

49.701 

(49.526) 

49.695 

(49.523) 

49.695 

(49.523) 

49.695 

(49.523) 

49.695 

49.695 

3.20 

11.266 

(10.873) 

11.265 

(10.872) 

11.265 

(10.372) 

11.265 

(10.872) 

11.265 

11.265 

6.40 

1.9774 

(1.6669) 

1.9767 

(1.6669) 

1.9767 

(1.6669) 

1.9767 

(1.6669) 

1.9767 

1.9767 

* 

CST  solution 


r 


r 
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Next,  the  results  of  bending  of  cross-ply  laminated  shells  are  discussed 
The  same  geometric  parameters  as  used  in  the  isotropic  shells  are  also  used 
for  composite  shells.  Individual  layers  are  assumed  to  be  orthotropic,  with 
the  following  properties: 

E^  =  25E2»  Gpj  =  0*2E,,»  6^3  =  Gj^  =  G.SE^,  =  D.25  (38) 

Tables  3-5  contain  nondimensional ized  center  transverse  deflections  of  cross- 
ply  spherical  shells  under  sinusoidal,  uniformly  distributed,  and  point  loads 
respectively.  The  results  are  tabulated  for  various  ratios  of  radius  to 
thickness  and  for  two  values  of  side  to  thickness.  From  these  results  it  fol 
lows  that  the  center  deflection  varies  rapidly  with  the  ratio  R/h  for  deep 
shells  (i.e.,  for  large  ratios  of  a/h)  than  for  shallow  shells  (i.e.,  for 
small  ratios  of  a/h). 

The  present  solution  can  also  be  applied  to  special  laminates, 
[0°/±45o]syrn  =  [0°/45o/-45o/-45o/45o/0°],  which  are  used  in  F-16  aircraft  by 
General  Dynamics  Corporation,  Ft.  Worth.  For  the  scheme,  ro°/±45°lsypi,  we 
have  Ajg  =  A2g  =  Ojg  =  =  8jj  =  f)  and  therefore  the  solution  developed  in 

the  present  study  applies.  Table  6  contains  the  nondimensional i zed  center  de 
flections  for  cylindrical  and  spherical  shells  with  various  side  to  thickness 
ratios  and  radius  to  side  ratios. 

To  complete  the  analysis,  the  results  of  free  vibration  are  also  present 
ed.  Tables  7  and  8  contain  nondimensional i zed  fundamental  frequencies  for 
spherical  (R]  =  R2'  and  ^ indrical  shells  (R]  =  1030)f  respectively. 

Two  sets  of  shear  correction  factors  were  used  to  investigate  their  in¬ 
fluence  on  the  fundamental  frequencies.  For  a  fixed  ratio  of  R/a,  the  shear 
correction  factors  have  little  or  no  effect  on  the  frequencies  for  a/h  =  100. 


V. 
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For  a/h  =  10,  the  effect  of  smaller  values  of  shear  correction  factors  is  to 
lower  the  frequencies,  the  rate  of  decrease  being  the  highest  for  the  four- 
layer  cross-ply.  For  identical  lamination  and  geometry,  spherical  shells  have 
higher  natural  frequencies  than  cylindrical  shells. 


★ 

TABLE  5.  Nondimensional ized  center  deflection  versus  radius  to  side  ratio  of 
spherical  shells  under  point  load  at  the  center  (101-term  solution) 


R/a 

0d/90° 

0<790!j/04  0  ^  90  a  7  9  0  0 

7o^"~ 

10 

4.0164 

3.7127  3.5815 

20 

5.7620 

4.6502  a. 4678 

30 

6.5444 

4.9546  4.7579 

40 

6.9098 

5.0794  4.8771 

50 

7.1015 

5.1410  4.9360 

100 

7.3836 

5.2273  5.0186 

1030 

7.4853 

5.2572  5.0472 

wh  JF„ 

w  =  (~ 

— )  x  102  ;  a/b  =  1 

,  a/h  =  10  ,  P  =  point  load 

Pa 

2 

TABLE  6. 

Nondimensional  center 

deflections  of  simply  supported 

,  special 

laminates  C0°/+45°]syrn 

cylindrical  and  spherical  shells. 

a/h 

Sinusoidal  loading  Uniform  loading 

* 

P/a 

5  10 

100  5  10 

ion 

Cylindrical  Shell  ( R / a  =  ?) 

2 

9.4976  4.1168 

0.1623  14.0224  6.1794 

0.0928 

a 

10.5152  4.9927 

0.5697  15.6616  7.5956 

0.733 

Spherical  Shell  =  R,  R/a  =  2) 

2 

6.7850  2.4120 

0.0421  9.7344  3.4794 

0.  0391 

4 

9.468  4.1117 

0.1623  13.9936  6.1856 

0.1930 

★ 


51-term  solution 
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TABLE  7  Nond i mens i onal i zed  fundamental  f requenci es ,+  w  =  wa 2/p/E,,/h ,  versus 
raaius  to  side  length  ratio  of  spherical  shell  (a/b=l). 


0° 

/90° 

b 

cn 

o 

o 

*/o° 

0  7  90° 

50 

O 

o 

O 

o 

R/a 

a/h=100 

a/h=10 

a/h=100 

a/h=10 

a/h= 100 

a/h=10 

1 

125.93 

14.481 

125.99 

16.115 

126.33 

16.172 

125.93 

14.420 

125.99 

15.746 

126.33 

15.827 

2 

67.362 

10.749 

68.075 

13.382 

68.294 

13.447 

67.361 

9.5161 

68.072 

12.794 

68.291 

12.8r0 

3 

46.002 

9.9608 

47.265 

12.731 

47.415 

12.795 

46.001 

9.6473 

47.260 

12.077 

47.410 

12.178 

4 

35.228 

9.4102 

36.971 

12.487 

37.082 

12.552 

35.227 

9.2644 

36.964 

11.808 

37.075 

11.910 

5 

28.825 

9.2309 

30.993 

12.372 

31.079 

12.437 

28.825 

9.0791 

30.986 

11.680 

31.071 

11.783 

10 

16.706 

8.9841 

20.347 

12.215 

20.380 

12.280 

16.704 

8.8225 

20.335 

11.506 

20.368 

11.609 

1030 

9.6873 

8.8998 

15.183 

12.162 

15.184 

12.226 

9.6850 

8.7317 

15.167 

11.447 

15.168 

11.551 

+the  first  line  of  values  corresponds  to  K?  =  K|  =  5/5,  and  the  second  line 
corresponds  to  K2  =  0.7,  K2  =  0.6, 


Table  8  Nondimensional ized  fundamental  frequencies  versus  radius  to  side 
length  of  cylindrical  shells  (a/b  =  1  ,  R-(  =  1 0  ^ u ) . 


0° 

/90° 

0°/90° 

/0° 

0°/90° 

/QQ°/0& 

Rp/a 

a/h=100 

a/h=10 

a/h=100 

a/h=10 

a/h=100 

a/h=10 

1 

65.474 

9.9986 

66.583 

13.172 

66.704 

13.128 

65.473 

9.9018 

66.580 

12.578 

66.700 

12.567 

;  2 

i 

34.914 

9.1476 

36.770 

12.438 

36.858 

12.471 

34.913 

9.0098 

36.763 

11.758 

36.851 

11.828 

|  3 

24.516 

8.9832 

27.116 

12.287 

27.173 

12.337 

24.515 

8.8832 

27.106 

11.587 

27.164 

11.676 

19.509 

8.9301 

22.709 

12.233 

22.749 

12.28Q 

19.508 

8.7746 

22.698 

11.526 

22.738 

11.622 

16.668 

8.9082 

20.332 

12.207 

20.361 

12.267 

b 

. 

16.667 

8.7498 

20.320 

11.498 

20.349 

11.596 

10 

11.831 

8.8879 

16.625 

12.173 

16.634 

12.236 

11.829 

8.7241 

16.610 

11.459 

16.619 

11.562 

1030 

9.6873 

8.8998 

15.183 

12.162 

15.184 

12.226 

(plate) 

9.6850 

8.7317 

15.167 

11.447 

15.168 

11.551 

t 
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CONCLUSIONS 

Closed-form  solutions  for  deflections  and  natural  frequencies  of  simply 
supported,  cross-ply  laminated  and  quasi -isotropic  shells  are  derived  using 
the  shear  deformation  version,  developed  herein,  of  the  Sanders  shell  theory. 
Unlike  plates,  antisymmetric  angle-ply  shells  do  not  admit  exact  solutions. 

The  exact  solutions  presented  herin  for  cylindrical  and  spherical  cross-ply 
shells  under  sinusoidal,  uniformly  distributed,  and  point  loads  should  serve 
as  bench  mark  results  for  approximate  methods,  such  as  the  finite  element  and 
finite  difference  methods. 
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APPENDIX  II.  -  NOTATION 

The  following  symbols  are  used  in  this  paper: 


d Aj »dA2 


F-i,E2 


=  extensional  stiffness;  see  Eq.  (21) 

=  curvilinear  dimension  of  the  laminate  along  the  £j-axis 
=  area  elements  perpendicular  to  the  ^  and  axes 
=  bending-extensional  stiffness;  see  Eq.  (21) 

=  curvilinear  dimension  of  the  laminate  along  the  ^?-axis 
=  coefficients  defined  in  Eq.  (29) 

=  plane-stress-reduced  material  stiffness  of  the  lamina 
=  bending  stiffnesses;  see  Eq.  (21) 

=  Young's  moduli  in  1  and  2  material  principal  directions 
=  shear  moduli  in  the  i-j  surfaces,  respectively 
=  thickness  of  the  shell 


=  Lame'  coefficients  (i  =  1,2,3) 

=  moment  resultants  (i  =  1,2,6);  see  Eqn.  (9) 

=  number  of  terms  in  the  series  (24)  and  (25) 

=  total  number  of  layers  in  the  laminate 
=  normal  stress  resultants  (i  =  1,2,6);  see  Eqn.  (9) 

=  number  of  terms  in  the  series  (24)  and  (26) 

=  rotary,  and  coupled-rotary  inertias  defined  in  Eqn.  (31) 

=  shear  stress  resultants  (i  =  1,2);  see  Eqn.  (9) 

=  coefficients  in  the  expansion  of  the  transverse  load;  see  Eq. 


=  position  vector  of  a  point  in  the  shell 
=  principal  radii  of  curvature;  see  Fig.  1 
=  position  vector  of  a  point  on  the  shell  midsurface 
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=  the  mass  coefficients  defined  in  Eqn.  (30) 

=  amplitudes  of  displacement  Uj ;  see  Eqn.  (25) 

=  displacements  of  the  midsurface 
=  displacements  of  a  point  in  the  shell 
=  amplitudes  of  displacement  u?;  see  Eqn.  (25) 

»  amplitudes  of  displacement  u^ 

=  amplitudes  of  rotation 

1  • 

=  amplitudes  of  rotation 
=  mu/a 

a  surface  metrics  (i  =  1,2)  defined  in  Eq.  (2) 

=  nu/b 

a  variational  operator 

=  rotations  about  normal  to  the  shell  midsurface 
=  strain  components  (i  =  1,2, 4, 5, 6);  see  Eqn.  (11) 

=  stress  components  (i  =  1,2, 4, 5, 6) 

=  curvilinear  coordinates  in  the  surface  of  the  shell 
=  coordinate  transverse  to  the  shell  midsurface 
=  Poisson's  ratios 
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ABSTRACT 

The  paper  contains  a  description  of  a  doubly  curved  shell  finite  element 
for  geometrical ly  nonlinear  (in  the  von  Karman  sense)  analysis  of  laminated 
(doubly-curved)  composite  shells.  The  element  is  based  on  an  extension  of  the 
Sanders  shell  theory  and  accounts  for  the  von  Karman  strains  and  transverse 
shear  strains.  The  numerical  accuracy  and  convergence  characteristics  of  the 
element  are  further  evaluated  by  comparing  the  present  results  for  the  bending 
of  isotropic  and  orthotropic  plates  and  shells  with  those  available  in  the 
literature.  The  many  numerical  results  presented  here  for  the  geomertical ly 
nonlinear  analysis  of  laminated  composite  shells  should  serve  as  reference  for 
future  investigations. 

INTRODUCTION 

Laminated  shells  are  finding  increased  application  in  aerospace,  automo¬ 
bile  and  petrochemical  industries.  This  is  primarily  due  to  the  high  stiff¬ 
ness  to  weight  ratio,  high  strength  to  weight  ratio,  and  less  machining  and 
maintenance  costs  associated  with  composite  structures.  However,  the  analysis 
of  composite  structures  is  more  complicated  when  compared  to  metallic  struc¬ 
tures,  because  laminated  composite  structures  are  anisotropic  and  character¬ 
ized  by  bending-stretching  coupling.  Further,  the  classical  shell  theories, 
which  are  based  on  the  Ki rchoff-Love  kinematic  hypothesis  (see  Naghdi  f 1 1  and 


Bert  [2])  are  known  to  yield  deflections  and  stresses  in  laminated  shells  that 
are  as  much  as  30%  in  error.  This  error  is  due  to  the  neglect  of  transverse 
shear  strains  in  the  classical  shell  theories. 

Refinements  of  the  classical  shell  theories  ( e. g. ,  Love's  first  approxi¬ 
mation  theory  [3])  for  shells  to  include  transverse  shear  deformation  have 
been  presented  by  Reissner  [4-6].  Sanders  [7]  presented  modified  first-  and 
second-approximation  theories  that  removed  an  inconsistency  (nonvanishing  of  a 
small  rigid-body  rotations  of  the  shell)  existed  in  Love's  first-approximation 
theory. 

The  first  thin  shell  theory  of  laminated  orthotropic  composite  shells  is 
due  to  Ambartsumyan  [fl,9].  In  these  works  Ambartsumyan  assumed  that  the  indi¬ 
vidual  orthotropic  layers  were  oriented  such  that  the  principal  axes  of  mate¬ 
rial  symmetry  coincided  with  the  principal  coordinates  of  the  shell  reference 
surface.  Dong,  Pister,  and  Taylor  [10]  presented  an  extension  of  Donnell's 
shallow  shell  theory  [11]  to  thin  laminated  shells.  Using  the  asymptotic  in¬ 
tegration  of  the  elasticity  equations,  Widera  and  Chung  [12]  derived  a  first- 
approximation  theory  for  the  unsymmetric  deformation  of  nonhomogeneous ,  aniso¬ 
tropic,  cylindrical  shells.  This  theory,  when  specialized  to  isotropic  mate¬ 
rials,  reduces  to  Donnell's  shell  theory. 

The  effects  of  transverse  shear  deformation  and  thermal  expansion  through 
the  shell  thickness  were  considered  by  Zukas  and  Vinson  [13].  Dong  and  Tso 
[14]  constructed  a  laminated  orthotropic  shell  theory  that  includes  transverse 
shear  deformation.  This  theory  can  be  regarded  as  an  extension  of  Love's 
first-approximation  theory  [3]  for  homogeneous  isotropic  shells.  Other  re¬ 
fined  theories,  specialized  to  anisotropic  cylindrical  shells,  were  presented 
by  Whitney  and  Sun  [15],  and  Widera  and  Logan  [16,17], 
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The  finite-element  analysis  of  layered  anisotropic  shells,  all  of  which 
are  concerned  with  bending,  stability,  or  vibration  of  shells,  can  be  found  in 
the  works  of  Schmit  and  Monforton  [18],  Panda  and  Natarajan  [19],  Shivakumar 
and  Krishna  Murty  [20],  Rao  [21],  Siede  and  Chang  [22],  Hsu,  Reddy,  and  Bert 
[23],  Reddy  [24],  and  Venkatesh  and  Rao  [25,26],  Recently,  Reddy  [27]  extend¬ 
ed  the  Sanders  theory  to  account  for  the  transverse  shear  strains,  and  pre¬ 
sented  exact  solutions  for  simply  supported  cross-ply  laminated  shells.  All 
of  these  studies  are  limited  small  displacement  theories  and  static  analyses. 

In  the  present  paper,  an  extension  of  the  Sanders  shell  theory  that  ac¬ 
counts  for  the  shear  deformation  and  the  von  Karman  strains  in  laminated  an¬ 
isotropic  shells  is  used  to  develop  a  displacement  finite  element  model  for 
the  bending  analysis  of  laminated  composite  shells.  The  accuracy  of  the  ele¬ 
ment  is  evaluated  by  comparing  the  results  obtained  in  the  present  study  for 
isotropic  and  orthotropic  plate  and  shell  problems  with  those  available  in  the 
literature.  Numerical  results  for  bending  analysis  of  cylindrical  and  doubly- 
curved  shells  are  presented,  showing  the  effect  of  radius-to-thickness  ratio, 
loading,  and  boundary  conditions  on  the  deflections  and  stresses. 

A  REVIEW  OF  THE  GOVERNING  EQUATIONS 
Consider  a  laminated  shell  constructed  of  a  finite  number  of  uniform¬ 
thickness  orthotropic  layers,  oriented  arbitrarily  with  respect  to  the  shell 
coordinates  The  orthogonal  curvilinear  coordinate  system 

(Sj.^.C)  is  chosen  such  that  the  and  ^ -  curves  are  lines  of  curvature  on 
the  midsurface  C=0,  and  C-curves  are  straight  lines  perpendicular  to  the  sur¬ 
face  C=0  (see  Fig.  1).  A  line  element  of  the  shell  is  given  by  (see  Reddy 
[27]_ 

(dS)2  *  [(1  +  c/R^ajdCj]2  +  [(1  +  m2)^2}2  *  (dC) 2  (1) 


'  -s— • 
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where  ai  and  R-j  (i  *  1,2)  are  the  surface  metrics  and  radii  of  curvature, 
respectively.  In  general,  and  R^  are  functions  of  only.  For  the  doubly 
curved  shells  considered  in  the  present  theory,  and  R,  are  constant. 

The  strain-displacement  equations  of  the  shear  deformable  theory  of 


doubly-curved  shells  are  given  by 


el  =  el  +  C<1 
-  o  ,  _ 

e«  -  €o  + 


e4  =  e4 


e5  =  e5 
o 

=  e6  +  C<6 


where 


L+Ol  +  L  (^1)2  K  s 

d  2  '■ax  *  i 


dUn 

-0  -  _ £  + 


l  *  H  .  i  &2 


Ki  =  ax. 


e6  =  dx. 


du^  3u^ 
+  9x1  ax^ 


a  $2  au2  aul 

k6  ”  ax,,  +  ax,  +  co  ^ ax,  ’  axT^ 


0  au3  u2 

e4  ‘  *2  +  ax^  "  Rj 


►l  +  W, 


r  -  1  r1  1  1 
C0  '  1  ^  “  R^> 


Here  u-j  denotes  the  displacements  of  the  reference  surface  along  =  C) 

axes,  and  are  the  rotations  of  the  transverse  normals  to  the  reference 
surface.  In  Love's  fi rst-approximation  theories  the  parameter  c0  is  taken  to 
be  zero,  and  it  is  introduced  only  in  the  Sanders  theory. 


,  t.  •  ki.* ’*■*>, . 
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The  stress-strain  relations,  transformed  to  the  shell  coordinates,  are  of 
the  form 

{a}  3  [Q] {e}  (4) 


(k) 

where  Qlj  are  the  material  properties  of  k-th  layer. 

The  principle  of  virtual  work  for  the  present  problem  is  given  by 


n  M  V  rCk  r  i  (k).  .  (kh  .  (k) _  .  (k)*  .  (k). 

0=2  J  /  loj  5ej  ♦  c£  6e2  +  6e6  +  '6e4  +  a £  '6 t 


k=l  C|M  Q 


q6Uj 


(5) 


=  /  +  N^Se^  +  +  M^6tc^  +  M26ic2  +  Mg6*g  +  Qj6eg 

+  Q26e^  -  q6u3]a1a2d^d52  (6) 


where  q  is  the  distributed  transverse  load,  and  are  the  stress  and 
moment  resultants,  and  0^  is  the  shear  force  resultant: 

L  ’’k 

(Ni,Mi)=  l  /  <*,  (l,C)dC  ,  1  «  1.2,6 

k=I  ck-l 
L  Ck 

Qi  »  I  Kf  /  c,dC  ,  i  =  4,5,  (7) 

ksl  ck=l 

where  (i  *  1,2)  are  the  shear  correction  factors  (taken  to  be  = 

5/6),  and  (Ck_i,Ck)  are  the  C -coordinates  of  the  k-th  layer,  and  L  is  the 
total  number  of  layers  In  the  laminated  shell. 

It  is  informative  to  note  that  the  equations  of  equilibrium  can  be 
derived  from  Eq.  (6)  by  integrating  the  displacement  gradients  in  e?  by  parts 
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and  setting  the  coefficients  of  6u.  (i  =  1,2,3)  and  6$.  (i=l,2)  to  zero 
separately.  We  obtain  [with  cQ  =  ^  {^—  -  -jj-)  and  dxi  =  cc^d!^  ] 

3N,  * 

4 ■ 0 

(N,  +  CM,)  +  +  jp  =  0 

3x^'  6  oo  3X2  ^2 


30 1  3Qp  Ni  Np 

hr*/-  q)  *  (“,)■« 


ax1  3x2 


«1  R2 


3M,  3M/r 

3^  +  3^  -  Q1  ’  ° 


3M,  3M„ 

3*f  +  8*f  -  *2  - 0 


(8) 


where 


-  3U-, 

C“n)  ■  af-  (Ni  »r*  N 


3U- 

'6  3x, 


3x, 


(N 


Su¬ 
fi  3x7 


3u- 


+  N?«f) 


(9) 


The  resultants  (Nit  Mi ,  0-j ) 


are  rel ated  to  ( e? ,  ) 


(i  .j 


1,2,6)  by 


Ni 

=  A..e° 
1J  J 

+ 

Bij'Cj 

Mi 

=  B .  .  e°- 
ij  J 

+ 

Vj 

°2 

=  A44eJ 

+ 

A45'5 

°1 

=  A45  e4 

+ 

A55  e5 

(10) 

(ID 


Here  A^j,  B^-  and  D^j  (i,j  =  1,2,6)  denote  the  extensional,  flexural  - 
extensional  coupling,  and  flexural  stiffnesses  of  the  laminate: 


(WV 


L  Ck  flcl 

y  J  0  1  (l,c,c2)dc  ( i ,j  =  1,2,6) 

k=1  ck-l 


(12) 
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(A44’A45’A55' 


L 

y 

k=l 


•k-1 


(*?>« 


•K1IC2Q45)  *K’°«' 


,00 

'55 


The  boundary  conditions,  derived  from  the  virtual  work  statement,  involve 
specifying  either  the  essential  boundary  conditions  (EBC)  or  the  natural 
boundary  conditions  (NBC): 


EBC 

U1 

u2 


u3 


♦l 


or 

or 

or 

or 

or 


■< 


me 

Vl  *  '"6  -  coM6>n2 

N?n2  *  (N6  *  coHfi)n2 

5Uo  dtio 

(N1  &x^nl  +  (N2  5x^n2 

du-  du., 

+  (N6  lix^nl  +  !'H6  ax^n2 


+  Q2n2  +  Qlnl 


Mlnl  +  M6n2 


M2n2  +  M6ni 


(13) 


where  (n^n^)  denote  the  direction  cosines  of  the  unit  normal  on  the  boundary 
of  the  midsurface  of  the  shell. 

The  exact  form  of  the  spatial  variation  of  the  solution  of  Eqs.  (8)-(13), 
for  the  smal 1 -di spl acement  theory,  can  be  obtained  under  the  following  condi¬ 
tions  (see  Reddy  [27]): 

(i)  Symmetric  or  antisymmetric  cross-ply  laminates:  i.e.,  laminates 
with 

a16  =  A26  =  B16  =  B26  =  °16  =  °26  =  A45  =  °* 

(ii)  Freely  supported  boundary  conditions: 

Nx(0,x2)  =  N1(a,x2)  =  M1(0,x2)  =  Mj(a,x2)  =  0 
u3(0,x2)  =  u3(a,x2)  =  u2(0,x2)  =  u2(a,x2)  =  0 


_7 
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N2(x1.0)  =  N2(xrb)  =  M2(xrO)  =  M2(xt,b)  =  0 

u3(xlf0)  =  u3(x1,b)  =  Uj(xpO)  =  u^Xj.b)  =  0 

4>2^ 0 »x 2 ^  ”  ^(a »x2^  ~  =  ^1^*1*^  “  ^  (15) 

where  a  and  b  are  the  dimensions  of  the  shell  middle  surface  along 

the  xj  and  X2  axes,  respectively.  The  time  variation  of  the  load 
does  not  influence  the  spatial  form  of  the  solution. 

Note  that  the  exact  solution  can  be  obtained  only  for  cross-ply  laminated 
shells  with  simply  supported  boundary  conditions.  For  general  lamination 
schemes,  exact  solutions  are  not  available  to  date. 

FINITE-ELEMENT  MODEL 

A  typical  finite  element  is  a  doubly-curved  shell  element  in  the 

/  Q  ) 

surface.  Over  the  typical  shell  element  Q  ,  the  displacements 
(uj,u2,u3,$j,$2)  are  interpolated  by  expressions  of  the  form, 

N  i 

u.  =  l  uf<b.(x  ,x  )  ,  i  =  1,8,3 

1  j=l  J  1  c 

N  i 

4us  I  ^(x,  ,x?)  ,  i  =  1,2  (16) 

1  j=l  J  1  i 

where  4,j  are  the  interpolation  functions,  and  u^  and  ^  are  the  nodal  values 
of  u-j  and  <t> • ,  respectively.  For  a  linear  isoparametric  element  (N  =  4)  this 
interpolation  results  in  a  stiffness  matrix  of  order  20  by  20.  For  a  nine- 
node  quadratic  element  t.ne  element  stiffness  matrix  is  of  order  45  by  45. 

Substitution  of  Eq.  (2)  into  the  virtual  work  principle,  Eq.  (9)  yields 
an  element  equation  of  the  fc rm 
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[K( A)  ]  =  {F}  (17) 

where  (a}  =  {{u^},  { } ,  {u^},  { ,  { 4)^ }  »  [K]  the  element  stiffness  matrix, 

and  { F }  is  the  force  vector.  In  the  interest  of  brevity,  the  coefficients  of 
the  stiffness  matrices  are  included  in  Appendix  I. 

The  element  equations  (17)  can  be  assembled,  boundary  conditions  can  be 
imposed,  and  the  resulting  equations  can  be  solved  at  each  load  step.  Note 
that  the  stiffness  matrix  [K]  is  a  function  of  the  unknown  solution  vector 
(a};  therefore,  an  iterative  solution  procedure  is  required  for  each  load 
step.  In  the  present  study,  we  used  the  direct  'teration  technique,  which  can 
be  expressed  as 

[K((slr)]{Alr+1  =  IF)  (18) 

where  (A}r  denotes  the  solution  vector  obtained  in  the  r-th  iteration  (at  any 
given  load  step).  At  the  beginning  of  the  first  load  step,  we  assume  that 
(A}°  =  {0}  and  obtain  the  linear  solution  at  the  end  of  the  first  iteration. 
The  solution  obtained  at  the  end  of  the  r-th  iteration  is  used  to  compute  the 
stiffness  matrix  for  the  (r+])-th  iteration.  At  the  end  of  each  iteration 
(for  any  load  step),  the  solutions  obtained  in  two  consecutive  iterations  are 
compared  to  see  if  they  are  close  enough  to  terminate  the  iteration  and  to 
move  on  to  the  next  load  step.  The  following  convergence  criterion  is  used  in 
the  present  study: 

1 1 1<  -  *r‘i2/ "  ki2iu?  <  3.01  (iq) 

i  =1  /  i  =  i 

where  N  is  the  total  number  of  unknown  generalized  displacements  in  the  finite 
element  mesh. 
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To  accelerate  the  convergence,  a  weighted  average  of  the  solution  from 
last  two  iterations  are  used  to  compute  the  stiffness  matrix: 

[K(,(A)r-1  *  (1  -  T){A)r)]Mrtl  =  |F|  (20) 

where  y  is  the  acceleration  parameter,  0  <  y  <  1.  In  the  present  study  a  val¬ 
ue  of  0.25  -  0.35  was  used. 


NUMERICAL  RESULTS 

Here  we  present  numerical  results  for  some  sample  problems.  To  illus¬ 
trate  the  accuracy  of  the  present  element,  first  few  examples  are  taken  from 
the  literature  on  isotropic  and  orthotropic  shells.  Then  results  ( i . . ,  de¬ 
flections  and  strsses)  for  several  laminated  shell  problems  are  presented. 

The  results  for  laminated  shells  should  serve  as  references  for  future  inves¬ 
tigations. 

All  of  the  results  reported  here  were  obtained  using  the  double-precision 
arithmatic  on  an  IBM  3081  processor.  Most  of  the  sample  problems  were  an¬ 
alyzed  using  a  2  x  2  uniform  mesh  of  the  nine-node  (quadratic)  isoparametric 
rectangular  element. 

1.  Bending  of  a  simply  supported  plate  strip  (or,  equivalently,  a  beam)  under 
uniformly  distributed  load. 

The  problem  is  mathematical ly  one-dimensional  and  an  analytical  solution 
of  the  problem,  based  on  the  classical  theory,  can  be  found  in  Timoshenko  and 
Womowsky-Krieger  [28],  The  plate  length  along  the  y-coordinate  is  assumed  to 
be  large  compared  to  the  width,  and  it  is  simply  supported  on  edges  parallel 
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to  the  y-axis.  The  following  simply  supported  boundary  conditions  are 
assumed : 

w  =  <jjy  =  0  along  edges  x  =  0,a.  (21) 

All  inplane  displacement  degrees  of  freedom  are  restrained.  A  5  x  1  mesh  of 
four-node  rectangular  elements  in  the  half  plate  is  used  to  analyze  the  prob¬ 
lem.  The  data  and  results  are  presented  in  Fig.  2.  The  present  result  is  in 
good  agreement  with  the  analytical  solution. 

2.  Clamped  square  plate  under  uniform  load. 

Due  to  the  biaxial  symmetry,  only  one  quadrant  of  the  plate  is  modelled 
with  the  2x2  mesh  of  nine-node  elements  (4x4  mesh  of  linear  elements  give 
almost  the  same  result).  Pertinent  data  and  results  are  presented  in  Fig.  3 
for  side  to  thickness  ratios  a/h  =  10  and  500.  The  result  for  a/h  =  500  is  in 
agreement  with  the  results  of  Way  [29],  The  difference  is  attributed  to  the 
fact  that  the  present  model  includes  the  inplane  displacement  degrees  of  free¬ 
dom  and  transverse  shear  deformation. 

Figure  4  contains  transverse  deflection  versus  load  for  clamped  ortho¬ 
tropic,  cross-ply,  and  angle-ply  plates.  The  lamina  properties  are 

E1  =  25  x  104  N/mm2,  E?  =  2  x  104  N/mm2,  G12  =  G13  =  104  N/mm2 

0^3  =  0.4  x  104  N/mm2,  =  0.25. 

For  the  same  total  thickness  the  clamped  orthotropic  square  plate  is  stiffer 
than  both  two-layer  angle-ply  and  cross-ply  plates. 
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3.  Simply  supported,  isotropic  spherical  shell  under  point  load. 

The  pertinent  data  of  the  shell  is  shown  in  Fig.  5.  A  uniform  mesh  of  2 
by  2  quadratic  elements  is  used  in  a  quadrant.  The  effect  of  three  types  of 
simply  supported  conditions  on  the  center  deflection  and  center  normal  stress 
is  investigated: 

SS-1:  u  =  w  =  4»=0  aty  =  b;  v  =  w  =  cj,=0  atx  =  a 

x  y 

SS-2:  u  =  v  =  w  =  4,  =  0  at  y  =  b;  u  =  v  =  w  =  cj>  =  0  atx  =  a  (22) 

x  y 

SS-3:  v  =  w  =  4j  =0  aty  =  b;  u  =  w  =  4,  =  0  at  x  =  a 

x  y 

Table  1  contains  the  results  for  the  three  boundary  conditions.  It  is  clear 
from  the  results  that  all  three  boundary  conditions  give  virtually  the  same 
results  for  a/h  =  160,  and  differ  significantly  (especially  SS-1  differs  from 
both  SS-2  and  SS-3)  for  a/h  =  16.  Thus,  the  effect  is  more  in  thick  shells 
than  in  thin  shells.  The  stress  Oy  shown  in  Fig.  5  is  evaluated  at  point  x  = 
y  =  1.691"  in  the  top  layer 

4.  Simply  supported  isotropic  cylindrical  shell  under  point  load. 

The  geometry  and  material  properties  of  the  shell  are  shown  in  Fig.  6. 

Once  again,  the  effect  of  various  simply  supported  boundary  conditions  (22)  on 
the  deflections  and  stresses  for  the  problem  is  investigated  using  a  uniform 
mesh  of  2  x  2  quadratic  eleme'  The  results  are  presented  in  Table  2.  For 
the  geometry  and  loading  used  here  (R  =  2540,  a  =  254,  h  =  12.7),  the  boundary 
conditions  have  very  significant  effect  on  the  solution.  Boundary  conditions 
SS-2  and  SS-3  give  almost  the  same  results  whereas  SS-1  gives  about  2  1/2 
times  the  deflection  given  by  SS-2  or  SS-3  boundary  conditions. 
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5.  Clamped  isotropic  cylindrical  shell  under  uniform  loading. 


Figure  7  contains  the  pertinent  data  and  results  for  a  clamped  cylindri¬ 
cal  shell  (isotropic)  subjected  to  uniform  load.  The  results  are  compared 
with  those  obtained  by  Ohatt  [30].  The  agreement  is  very  good. 

6.  Clamped  orthotropic  cylindrical  shell  subjected  to  internal  pressure. 

Figure  8  contains  the  geometry  and  plots  of  center  deflection  and  center 
stress  versus  the  internal  pressure  for  the  problem.  The  orthotropic  material 
properties  used  in  the  present  study  are: 

E1  =  7.5  x  106  psi ,  =  2  x  106  psi,  =  G13  =  =  1 .25  x  106  psi 

v12  =  0.25  (23) 

The  present  result,  obtained  using  the  2x2  mesh  of  quadratic  elements,  is  in 
excellent  agreement  with  that  obtained  by  Chang  and  Sawamiphakdi  [31]. 

7.  Nine-layer  [0o/90°/0°. . ,/0°]  cross-ply  spherical  shell  subjected  to 
uniformly  distributed  load. 

The  following  geometrical  data  is  used  in  the  analysis: 

Rl  =  =  R  =  1,000  in.,  a  =  b  =  100  in.,  h  =  1  in.  (24) 

Individual  layers  are  assumed  to  be  of  equal  thickness  (h^  =  h/9),  with  the 
zero-degree  layers  being  the  inner  and  outer  layers.  The  following  two  sets 
of  orthotropic-material  constants,  typical  high  modulus  graphite  epoxy  materi¬ 
al  (the  ratios  are  more  pertinent  here),  for  individual  layers  are  used: 
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Mat.-l:  E}  =  25  x  106  psi ,  E2  =  106  psi ,  G12  =  G13  =  0.5  x  106  psi 

G23  =  0.2  x  106  psi,  v12  =  0.25  (25) 

Mat. -2:  E^  =  40  x  106  psi,  E2  =  106  psi,  G^2  =  Gj3  =  0.6  x  106  psi 

G23  =  0.5  x  106  psi,  v12  =  0.25  (26) 

Figure  9  contains  plots  of  center  deflecton  (w/h)  versus  the  load  parameter 
(P  =  qQR2/E2h2)  for  the  two  materials.  Shell  constructed  of  Material  1 
deflects  more,  for  a  given  load,  than  the  shell  laminated  of  Material  2 
(because  Material  2  is  stiffer),  and  consequently  experiences  greater  degree 
of  nonlinearity.  Note  that  the  difference  between  the  nonlinear  deflections 
of  the  two  shells  increase  nonlinearly,  indicating  that  the  shell  made  of 
Material  2  can  take  much  more  (ultimate)  load  than  apparent  from  the  ratio  of 
moduli  of  the  two  materials,  e[^/e|^. 

8.  Effect  of  various  simply-supported  boundary  conditions  on  the  deflections 
of  two-layer  cross-ply  spherical  shells  under  uniform  load. 

As  pointed  out  in  Problems  3  and  4,  the  transverse  deflection  is  sensi¬ 
tive  to  the  boundary  conditions  on  the  inplane  di spl acements  of  simply  sup¬ 
ported  shells.  To  further  illustrate  this  effect  for  laminated  shells,  a  set 
of  four  types  of  boundary  conditions  are  used,  and  the  results  are  presented 
in  Table  3.  Here  SS-4  has  the  following  meaning: 

w  =  i  =  0  on  x  =  0  and  a 

SS-4  (27) 

w  =  ^  =  0  on  y  =  0  and  b 
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Once  again  we  note  that  SS-2  and  SS-3  give  almost  the  same  deflections. 
Boundary  conditions  SS-1  and  SS-4  give  deflections  an  order  of  magnitude  high¬ 
er  than  those  given  by  SS-2  and  SS-3.  Thus,  boundary  conditions  SS-2  and  SS-3 
make  the  shell  quite  stiffer. 

9.  Two-layer  cross-ply  C0o/90o 3  and  angle-ply  [-45°/45°],  simply-supported 
(SS-3)  spherical  shells. 

Figure  10  contains  the  pertinent  data  and  results  (with  different  scales) 
for  the  cross-ply  and  angle-ply  shells  (of  Material  2).  It  is  interesting  to 
note  that  the  type  of  nonlinearity  exhibited  by  the  two  shells  is  quite  dif¬ 
ferent;  the  cross-ply  shell  gets  softer  whereas  the  angle-ply  shell  gets 
stiffer  with  an  increase  in  the  applied  load.  While  both  shells  have  bending¬ 
stretching  coupling  due  to  the  lamination  scheme  (822  =  "  Bn  nonzero  for  the 
cross-ply  shell  and  8^5  and  B2§  are  nonzero  for  the  angle-ply  shell),  the 
angle-ply  experiences  shear  coupling  that  stiffens  the  spherical  shell  rela¬ 
tively  more  than  the  normal  coupling  (note  that,  in  general,  shells  get  softer 
under  externally  applied  inward  load). 

Figure  11  contains  plots  of  center  deflection,  normal  stress  (-ay)  and 
shear  stress  (oy2)  at  x  =  y  =  5.283"  versus  load  for  two-layer  cross-ply 
(0°/90°)  spherical  shell  (Material  1)  under  point  load  at  the  center  of  the 
shell.  The  nonlinearity  exhibited  by  the  stresses  (especially  a  )  is  less 
compared  to  that  exhibited  by  the  transverse  deflection. 

10.  Two-layer  clamped  cylindrical  shells  under  uniform  loads. 

Figures  12  and  13  contain  results  (i.e.,  w ,  a  ,  a  versus  load)  for 

y  x  ^ 

cross-ply  [0°/90°]  and  angle-ply  [-45°/45°]  clamped  cylindrical  shells  under 
uniform  load.  The  load-deflection  curve  for  the  cros-ply  shell  resembles  that 
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of  the  isotropic  shell  in  Fig.  7,  but  exhibits  greater  degree  of  nonlinearity 
(being  stiffen).  The  angle-ply  shell  exhibits  different  type  of  nonlinearity 
(softening  type)  for  all  loads. 


11.  Quasi -i sotropic ,  clamped,  cylindrical  shell  under  uniform  load. 

Two  types  of  quasi-isotropic  clamped  cylindrical  shells  are  analyzed: 

Type  1:  [0°/45790°/-45°] 

y  •  (28) 

Type  2:  [0°/±45o/90]sym# 

Material  1  properties  are  assumed  for  each  lamina  (8  layers).  The  geometric 
data  and  results  are  presented  in  Fig.  14.  Compared  to  the  results  presented 
iri  Figs.  12  and  13,  the  quasi -i sotropic  shells  have  the  'near-inflection' 
point  at  higher  loads;  the  load -deflect ion  curve  has  essentially  the  same  form 
as  that  of  the  cross-ply  shell  (see  Fig.  12). 


SUMMARY  AND  CONCLUSIONS 

A  shear-flexible  finite  element  based  on  the  shear  deformation  version  of 
the  Sanders'  theory  and  the  von  Karman  strains  is  developed,  and  its  applica¬ 
tion  to  isotropic,  orthotropic,  and  laminated  (cross-ply  and  angle-ply)  shells 
is  illustrated  via  numerous  sample  problems.  Many  of  the  results,  especially 
those  of  laminated  shells,  are  not  available  in  the  literature  and  therefore 
should  serve  as  references  for  future  investigations.  From  the  numerical  com¬ 
putations  it  is  observed  that  boundary  conditions  on  the  inplane  displacements 
have  significant  effect  on  the  shell  deflections  and  stresses.  Also,  it  is 
noted  that  the  form  of  nonl inearities  exhibited  by  different  lamination 
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schemes  (i.e.,  cross-ply  and  angle-ply)  is  different.  It  would  be  of  consid¬ 
erable  interest  to  verify  these  findings  by  experiments. 
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Table  1.  Effect  of  various  simply  supported  boundary  conditions  on  the  center 
deflections  and  normal  stress  in  spherical  shells  under  point  load 
(E  =  107  psi,  v  =  0.3). 


Load  Solution  _ SS-1 _  S5-Z  rrlr - - 

p/h 2  a/h= 160  a/h*16  a/h=16(5 a/h=l6  a/h=160  a/h=16 


4,000 

-w* 

0.0155 

- 

-a  * 

893 

- 

8,000 

-wX 

0.0329 

0.0349 

-oy 

1,880 

6,535 

12,000 

X 

-w 

0.0529 

- 

2,980 

- 

16,000 

V 

0.0760 

0.0793 

-crv 

4,220 

13,230 

20,000 

X 

-w 

0.1038 

- 

-a 

5,657 

- 

24,000 

X 

-w 

0.1364 

0.1083 

-ay 

7,268 

20,110 

28,000 

X 

-w 

0.1761 

- 

-cy 

9,128 

- 

32,000 

X 

-w 

0.2234 

0.1472 

-°x 

11,180 

27,170 

0.0152 

- 

0.0152 

- 

984 

- 

894 

- 

0.0324 

0.0255 

0.0324 

0.0258 

1,882 

6,015 

1,882 

6,031 

0.0522 

- 

0.0521 

- 

2,985 

_ 

2,986 

- 

0.0752 

0.0520 

0.0751 

0.0525 

4.228 

12,200 

4,229 

12,240 

0.1028 

- 

0.1027 

- 

5,671 

- 

5,672 

- 

0.1354 

0.0792 

0.1353 

0.0800 

7,289 

18,500 

7,291 

18,550 

0.1752 

- 

0.1751 

- 

9,160 

- 

9,162 

- 

0.2227 

0.1072 

0.2227 

0.1083 

11,220 

24,930 

11,230 

25,000 

*  w( 0,0) ,  ox( A, A);  A  =  1.691 


Table  2.  Effect  of  various  types  of  simply  supported  boundary  conditions  on 
the  deflections  and  stresses  of  anisotropic  cylindrical  shell  under 
point  load. 


Load  ,P 
(N) 

/  \  ss- 

-w(mm) 

l-0y(N/mm2) 

SS-2 

-w 

•°y 

-w  SS-3 

■°y 

250 

2.5804(2) 

2.868 

0.6544(4) 

1.706 

0.6698(4) 

1.706 

500 

5.1626(2) 

5.713 

1.3533(4] 

i 

3.478 

1.3843(4) 

3.477 

750 

7.7343(2) 

8.506 

2.1057(4 

■ 

5.327 

2.1522(4) 

5.321 

1,000 

10.278(2) 

11.210 

2.9234(4] 

1 

7.265 

2.9855(4) 

7.242 

1,250 

12.733(2) 

13.80 

3.8241(4) 

i 

9.312 

3.9017(4) 

9.288 

1,500 

15.204(2) 

16.25 

4.8349(4] 

) 

11.50 

4.9279(4) 

11.46 

1,750 

17.560(2) 

18.560 

6.0331(5] 

1 

13.91 

6.1423(5) 

13.85 

2,000 

19.843(2) 

20.730 

7.5316(6' 

) 

16.66 

7.6610(6) 

16.57 

r  - »  A  ..  yp.T.  - 
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Table  3.  Effect  of  various  simply-supported  boundary  conditions  on  the  trans 
verse  deflections  of  cross-ply  [0°/90°]  spherical  shells  under 
uniform  load. 


9o 

SS-1 

SS-2 

SS-3 

SS-4 

0.50 

0.3344 

0.04246 

0.04257 

0.4592 

0.75 

0.5757 

0.06599 

0.06617 

0.8255 

1.00 

0.9485 

0.09144 

0.09171 

1.3845 

1.25 

1.6529 

0.11926 

0.11966 

1.9589 

1.50 

2.2826 

0.15008 

0.15063 

2.3597 

1.75 

2.6421 

0.18478 

0.18556 

2.5951 

2.00 

2.8499 

0.22473 

0.22584 

2.8074 

2.25 

3.0764 

0.27425 

0.27593 

3.0284 

2.50 

3.2432 

0.33534 

0.33795 

3.1948 

2.75 

3.4214 

0.42970 

0.43487 

3.3719 
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Stiffness 

Let 

[Kl1]  =  A 

[Kl2]  =  a 

[KI3]  =  f 


[K14]  =  B 


APPENDIX  I 

Coefficients : 

.  _  1_  ^3  f  _  1  0u3 

1  2  dxj  *  2  "  2  dx2 

iCSll]  +  A16([S12]  +  [S21])  +  A66[S22] 

c0(B16([S12]  .  [S21])  t  2b66[s22)  .  co066[s«])  *  ^5  [S°°] 

R1 

2[s12]  +  A16[Sll]  +  A26[S22]  +  A66[S21] 
co^26^s32^  "  B16^SU^  *  coD66I-S21^  * 

jfAjjCSH]  *  «16[S12]  *  CS2'])  *  Vs22]) 
f2A12[S‘n  ♦  A16[S“]  *  A^S22]  t  A6sCS2‘1) 
i-(4ll[sl0]+A16[S2O])tl-(A12[S*»]*A26[S2»]) 
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[K3i]  =  [2K13]^  ,  NL  =  Nonlinear  portion  of  the  matrix 
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r  kra^l  =  rt(f?alT 

jlinear  LN  JLinear 
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i  j 


-  / 


dx 
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dx^dx^ 


.  Sfj  =  /  4>142jdx1dx2 

Qe 


It  should  be  noted  that  although  f1  and  f2  are  shown  factored  outside  the 
matrices,  in  the  evaluation  of  the  coefficients  by  the  Gauss  quadrature  f^ 
f2  are  considered  as  parts  of  the  integrals.  For  example  fiA^Su]  is 


evaluated  by 
Q6 


,  N  N  du-j 

j  l  l  _7  x  -  W  W  ,detJ 

2  i-i  j=i  11  5xi  1  J  xrzrx2~zj  1  J 


where  N  is  the  numtei  of  Gauss  points,  Wj  and  Wj  are  the  Gauss  weights,  Zj 
Zj  are  the  Gauss  points,  and  J0  is  the  Jacobian  of  the  transformation. 
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Figure  3.  Bending  of  clamped  isotropic  square  plate  under 
uniform  load. 


4  « 


31 


Figure  4. 


[•TT« 


34 


Figure  7. 


Bending  of  a  clamped,  isotropic,  cylindrical  shell  under 
uniform  load. 
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Deflection,  (-w/h) 


Figure  9.  Bending  of  nine-layer  cross-ply 
[07  90  7  08/...]  spherical  shell 
subjected  to  uniformly  distributed 
load. 
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Figure  10.  Bending  of  two-layer  cross-ply  and  angle-ply, 
sinply  supported  (SS-3)  spherical  shells  under 
uniform  load. 
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Figure  11.  Bending  of  a  cross-ply  [0°/90°]  spherical  shell 
(SS-3,  Material  1),  under  point  load. 
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Figure  13. 


Bending  of  a  damped  angle-ply  [-450/45°l 
shell  under  uniform  load. 


cylindrical 


Figure  14.  Bending  of  clamped  quasi-isotropic  cylindrical  shells 
under  uniform  load. 


r 

UNCLASSIFIED 


SECURtTY  CLASSIFICATION  OP  THIS  PAGE  (When  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

t.  REPORT  NUMBER 

VPI-E-83. 45 

2.  GOVT  ACCESSION  NO. 

fib  -  fins 8<f? 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  T|Tl_E  (and  Subtitle) 

EXACT  AND  FINITE-ELEMENT  ANALYSIS 

SHELLS 

3F  LAMINATED 

S.  TYPE  OF  REPORT  4  PERIOD  COVERED 

Interim 

S.  PERFORMING  ORG.  report  NUMBER 

Tech.  Report  No.  AFOSR-81-5 

7.  autmo«c»; 

8.  contract  OR  grant  NUMBERfjj 

J.  N.  Reddy 

AF0SR-81  -01 42 -B  and  C 

9  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Virginia  Polytechnic  Institute  and  State  Universv 
Blacksburg,  Virginia  24061 

to.  program  ELEMENT.  PROJECT  task 
AREA  4  WORK  UNIT  NUMBERS 

2307/BI 

(Structural  Mechanics) 

n.  controlling  office  name  ano  aooress 

United  States  Air  Force 

Air  Force  Office  of  Scientific  Research 

Bldg.  410,  Boling  AFB,  D.C. 

12.  REPORT  OATE 

November,  1983 

13.  NUMBER  OF  PAGES 

70 

14.  monitoring  AGENCY  NAME  4  AOORESS (it  diftarent  trom  Controlling  Ollict) 

IS.  SECURITY  CLASS,  (of  this  report) 

Unclassified 

IS#.  OECL  ASSIPICATION.  OOwnGR  AOinG 

schedule 

16.  DISTRIBUTION  STATEMENT  (of  this  Report) 


This  document  has  been  approved  for  public  release  and  sale;  distribution 
uni  i mi  ted. 


17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  in  Black  20,  if  different  from  Report) 


is.  supplementary  notes 

The  report  has  two  parts.  The  first  part  deals  with  the  exact  solution  of 
the  linear  theory  governing  shear  deformable,  doubly-curved  laminated  shells. 
The  second  part  (coauthored  with  K.Chandrashekhara)  deals  with  the  finite-ele- 
— ment  analysis  of  the  same_ theory,  but  also  iadmies  the  wn  ^rman  _ 

19.  KEY  WOROS  (Continue  on  reverse  aid#  if  necessary  and  identify  by  block  number) 

Composite  materials,  exact  solutions,  finite  element  solutions,  laminates, 
plates,  shells,  shear  deformation  theory,  von  Kerman  strains  (geometric 
nonlinearity),  vibration  frequencies. 


20.  ABSTRACT  ( Continue  on  reverse  side  if  necessary  and  identify  by  block  number ) 

^The  two  parts  of  the  report  deal  with  the  exact  solutions  of  the  linear  theory 
and  finite-element  analysis  of  the  geometrically  nonlinear  (i.e.  the  von  Karman 
theory  of  laminated,  doubly-curved,  shells,  respectively.  Numerical  results 
are  presented  for  both  plates  and  shells.  The  results  should  serve  as  referent 
for  future  investigations. 

— — - ' 

DO  ,  j*nM73  1473  edition  a?  t  „ovm  .sossoLETe  UNCLASSIFIED 

CLASSIFICATION  OP  This  page  When~Datm  Entered' 


